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The concept of gauge invariance is one of the most subtle and useful concepts in modern theoretical 
physics. It is one of the Standard Model cornerstones. The main benefit due to the gauge invariance 
is that it can permit the comprehension of difficult systems in physics with an arbitrary choice of a 
reference frame at every instant of time. It is the objective of this work to show a path of obtaining 
gauge invariant theories from non-invariant ones. Both are named also as first- and second-class 
theories respectively, obeying Dirac's formalism. Namely, it is very important to understand why 
it is always desirable to have a bridge between gauge invariant and non-invariant theories. Once 
established, this kind of mapping between first-class (gauge invariant) and second-class systems, 
in Dirac's formalism can be considered as a sort of equivalence. This work describe this kind of 
equivalence obtaining a gauge invariant theory starting with a non-invariant one using the symplectic 
embedding formalism developed by some of us some years back. To illustrate the procedure it was 
analyzed both Abelian and non-Abelian theories. It was demonstrated that this method is more 
convenient than others. For example, it was shown exactly that this embedding method used here 
does not require any special modification to handle with non-Abelian systems. 

PACS numbers: 11.15.-q; ll.lO.Ef; 11.30.Cp 
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I. INTRODUCTION 

In this work the symplectic embedding of second-class 
systems transforming them into gauge theories is dis- 
cussed [1]. Although there is an extensive literature 
about the subject, we believe that the results, the gauge 
invariance and the symplectic embedding applications de- 
scribed here are shown in a self-contained and pedestrian 
way. 

The connection between these conversion features can 
be understood as a kind of physical equivalence. This 
connection characterizes both systems as representing 
the same physical properties @. The term embedding, 
used here (and in the literature about the subject) sum 
up precisely the gist of the procedure. We will use a 
general canonical formalism of embedding developed by 
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some of us based on symplectic formalism 3] , which em- 
beds a second-class theory (the initial action, noninvari- 
ant) into one that has gauge invariance. It is important 
to emphasize here that our embedding leads to a mixed 
system having both first and second class constraints. 

To perform such conversion task we then used the 
so-called symplectic embedding formalism developed by 
some of us in The main advantage of this method, 
when compared with other conversion methods, is that it 
involves the construction of a zero-mode object, which is 
connected with the symmetries of the final model. Hence, 
the embedding method discloses all the hidden symme- 
tries of the original system, since it can be demonstrated 
that the original and final models present the same equa- 
tions of motion 2\ . The revealing of the set of symmetries 
of the model is a fundamental issue concerning the fath- 
oming of the original action. Another advantage is that 
the method can deal with Abelian and non-Abelian the- 
ories. Namely, nothing needs to be modified or adapted 
when we use the embedding method to deal with such 
theories. In the last two sections we will show that this 
method can be used by a wider community since we at- 
tack fluid systems. 
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As well known, after the procedure application, where 
we have first-class constraints, the final system is a gauge 
invariant theory. The relevance of having a gauge the- 
ory, in few words, is that the physically important vari- 
ables are those that are independent of the local reference 
frame |4|. Whenever a change in the arbitrary reference 
frame causes a transformation of the variables involved 
we have the so-called gauge transformation. This physi- 
cal variables are then well known as gauge invariant vari- 
ables. Such gauge theories and gauge transformations 
are the cornerstones of the standard model structure, to 
mention only one of its successful applications. We will 
talk with more detail about gauge theories and quanti- 
zation in a moment. 

Since the Hamiltonian formulation is considered by 
many as the more fundamental formulation of a phys- 
ical theory, whenever necessary, we also used the Dirac 
brackets [l| in order to show exactly the gauge invariance 
of the final actions obtained. 

Quantization. The quantization of gauge theories 
demands a special care because the presence of gauge 
symmetries indicates that exist some superfluous degrees 
of freedom, which must be eliminated (before or after) 
with the implementation of a valid quantization process. 

The quantization of first-class systems was formulated 
both in Dirac's [l[ and path integral [5j points of view. 
Later on, the path integral analysis was extended by 
Batalin, Fradkin and Vilkovisky [6] in order to preserve 
the BRST symmetry Q. 

On the other hand, the covariant quantization of 
second-class systems is, in general, a difficult task be- 
cause the Poisson brackets are replaced by Dirac brack- 
ets. At the quantum level, the variables become opera- 
tors and the Dirac brackets become commutators. Hence, 
the canonical quantization process is contaminated with 
serious issues such as ordering operator problems and 
anomalies 0] in the context of nonlinear constrained sys- 
tems and chiral gauge theories, respectively. It seems 
that it is more natural and safe to work out the quan- 
tization of second-class systems without invoking Dirac 
brackets. Actually, it was the strategy followed by many 
authors over the last decades. The noninvariant system 
has been embedded in an extended phase space in or- 
der to change the second-class nature of constraints into 
first-class. 

In this way, the entire machinery 0, HH for quan- 
tizing first class systems can be used. To implement 
this concept, Faddeev [l4[ suggests to enlarge the phase 
space with the introduction of new variables to linearize 
the system, which were named after them, as the Wess- 
Zumino (WZ) variables [H)]. This idea has been em- 
braced by many authors and some methods were pro- 
posed and some constraint conversion formalisms, based 
on the Dirac method [l| , were constructed. Among them, 
the Batalin-Fradkin-Fradkina-Tyutin (BFFT) [10|, ED] 
and the iterative [l6| methods were strong enough to be 
successfully applied to a great number of important phys- 



ical systems. Although these techniques share the same 
conceptual basis [l4| and follow the Dirac framework [l[ , 
these constraint conversion methods were implemented 
following different directions. Historically, both BFFT 
and iterative methods were introduced to deal with lin- 
ear systems such as chiral gauge theories [HI, [ljj in order 
to eliminate the gauge anomaly that hampers the quan- 
tization process. 

The Symplectic formalism. The technique formu- 
lated by Faddeev and Jackiw [18| relies on first-order La- 
grangians which equations of motion does not generate 
accelerations. The brackets involved are called gener- 
alized brackets since it can be shown that they coin- 
cide with those obtained directly from Dirac's formal- 
ism. So, we can realize that these generalized brackets 
are linked to the commutators of the quantized theory. In 
FJ method the classification of a system as constrained 
or unconstrained is connected with the behavior of the 
symplectic two-form. The outcome of this classification 
does not necessarily agree with the one given by the Dirac 
technique. 

The BFFT technique. The general canonical quan- 
tization formalism due to Batalin, Fradkin, Fradkina 
and Tyutin (BFFT) [lOj for converting second-class con- 
straints into first-class ones, increases the number of de- 
grees of freedom to include unphysical ones 11]. The 
obtainment of first-class constraints is carried out in an 



iterative process 12[ . To summarize we can say that the 
first correction of the constraints is linear in the auxil- 
iary variables. The second correction is quadratic and 
so on. For systems with just linear constraints we have 
to perform the conversion to first-class ones. For nonlin- 
ear constraints, we have more than one iteration and the 
BFFT formalism shows its problems since the first step 
of the process does not precisely fix the solution that will 
be used in the next steps. Some solutions of this problem 
were published but they cannot be applied for all cases. 

The Noether dualization method. Recently, the 
so-called gauging iterative Noether dualization method 
[l9| has been shown to thrive in establishing some du- 
alities between models [2(|. This method hinges on the 
traditional concept of a local lifting of a global symme- 
try and it may be realized by an iterative embedding of 
Noether counterterms. However, this method provides a 
strong suggestion of duality since it has been shown to 
give the expected result in the paradigmatic duality be- 
tween the so-called self-dual model [21| and the Maxwell- 
Chern-Simons theory in three dimensions duality. This 
correspondence was first established by Deser and Jackiw 
[22| and using the parent action approach (23j. 

The paper. In section 2, we review the symplectic 
embedding formalism in order to settle the notation and 
to make the reader to familiarize with the fundamentals 
of the formalism. In section 3, we will initiate with the 
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Proca model in order to set up the general ideas dis- 
cussed in Section 2. We will see that a different zero- 
mode can bring us an equivalent action different from 
the literature. After that, we apply this formalism in 
two important models in high energy physics. The first 
one, in section 4, is the nonlinear sigma model (NLSM) 
[24| , which is an important theoretical laboratory to learn 
the basics about asymptotically free field theories, as dy- 
namical mass generation, confinement, and topological 
excitations, which is expected in the realistic world of 
four-dimensional non-Abelian gauge theories. The sec- 
ond, in section 5, is the bosonized chiral Schwinger model 
(CSM), which has attracted too much attention over the 
last decade, mainly in the context of string theories p9l |. 
and also due to the huge progress in the understanding of 
the physical meaning of anomalies in quantum field theo- 
ries achieved through the intensively study of this model. 
The gauge invariant version of the model was obtained. 
Section 6 is devoted to an application of the formalism 
to the non-Abelian Proca model. We will show that this 
gauge-invariant formalism does not require modifications 
to deal with non-Abelian models as demanded by the 
BFFT method. In section 7, we enlarged the method 
boundaries performing the embedding of a fluid dynami- 
cal model, showing that the method is not restricted only 
to the theoretical ones. In section 8, we analyze the sym- 
metries of the rotational fluid model with a new extra 
term, which introduced a dissipative force into the sys- 
tem. The objective is to promote an approximation to 
reality where we always have dissipation. Whenever con- 
venient, the gauge symmetry will be investigated from 
the Dirac point of view. Our concluding observations 
and final comments are given in Section 9. 



II. THE SYMPLECTIC EMBEDDING 
FORMALISM 

The formalism is developed based on the symplectic 
framework [l8l . [25| , that is a modern way to handle with 
constrained systems. The basic object behind this for- 
malism is the presymplectic matrix, i.e., if this matrix is 
singular, the model presents a symmetry, if it is not sin- 
gular, the Dirac brackets can be obtained. In this way, 
we have proposed to turn nonsingular presymplectic ma- 
trix to a singular one. This procedure will be carried out 
through the introduction of both the arbitrary functions 
that depend on the original coordinates and the WZ vari- 
ables into the first-order Lagrangian. To appreciate this 
point, a brief review of the symplectic formalism will be 
provided in this section. After that, general ideas of the 
symplectic gauge-invariant formalism will be explained. 
This formalism, differently from the BFFT and other it- 
erative constraint conversion methods, does not require 
special modifications of the formalism to convert Abelian 
or non-Abelian sets of second-class constraints into first- 
class ones. We think that this is the main advantage of 
this formalism. 



In the following lines, as we said above, we will try to 
keep this paper self-contained. We will follow closely the 
ideas contained in (26j . 

Let us now consider a general noninvariant mechani- 
cal model whose dynamics is governed by a Lagrangian 
hi, t), (with i = 1,2,..., N), where at and hi are the 
space and velocity variables, respectively. Following the 
symplectic method the zeroth-iterative first-order one- 
form Lagrangian is written as 



(2.1) 



with arbitrary one form AW = A ( g>d^ e and the 
presymplectic variables are defined as 



(0)0 



a l , with 6 = 1,2, N 

Pi, with e = N + 1, N + 2, . . . , 2N, 



(2.2) 



where A^ are the canonical momenta and V"' ' is defined 
as being the presymplectic potential. From the Euler- 
Lagrange equations of motion we have that 

dV 



(0)/(0)/3 = 



JOB S 



whose solutions rely on the invertibility of fgp. The 
problem is that it is an impossible task if there are con- 
straints involved. They would provide a singular matrix 
since 



(0)/3 



(f (oh-i dV 

> Q£(0)6 



(2.3) 



It is easy to see that the Hamiltonian corresponding to 
the Lagrangian (|2.ip is l/ (0) . So, Eqs. (|23|) are also 
Hamiltonian- type 



£(0)/3 _ {y(0)^(0)/3| 



dV 



{^(o)e^(o)/3i 



But the generalized bracket is defined to be 
£(0)/?} = (f^)- 1 . The equations of motion 
(12.3[) can be written in Hamiltonian form, with the 
presymplectic potential V^ ' playing the role of Hamil- 
ton function. 

The presymplectic tensor is obtained as 



Co) 



dA 



(o) 



dAf 



(2.4) 



If the two-form 



/ ee -f e ^ e A d£" 



is singular, the presymplectic matrix (|2.4[) has a zero- 
mode that generates a new constraint when mul- 
tiplied by (contracted with) the gradient of the presym- 
plectic potential, 



,(0)0 



(2.5) 
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This constraint will be introduced into the zeroth- with 
iterative one-form Lagrangian equation (|2.1[) through a 
Lagrange multiplier 77, generating the next one 

C^dt = A^d^ 9 + dtflW - V^iOdt, 

= A^d^-V^{Odt, (2.6) 

with 7 = l,2,...,(2iV + l) and 

ph = tfW, v)i (2.7) 

4> = (4 0) ,fi(°)). 



As a consequence, the hrst-iterative presymplectic tensor 
can be computed as 



dA 



(i) 



(2- 



If this tensor is nonsingular, the iterative process stops 
and the Dirac brackets among the phase space variables 
are obtained from the inverse matrix (/^) _1 and con- 
sequently, the Hamiltonian equation of motion can be 
formulated and solved, as discussed in [27]. It is well 
known that a physical system can be described at least 
classically in terms of a presymplectic manifold M. . From 
a physical point of view, M. is the phase space of the sys- 
tem while a nondegenerate closed 2-form / can be iden- 
tified as being the Poisson bracket. The dynamics of the 
system is determined just specifying a real- valued func- 
tion (Hamiltonian) H on the phase space, i.e., one of 
these real-valued function solves the Hamiltonian equa- 
tion, namely, 



i(X)f = dH, 



(2.9) 



and the classical dynamical trajectories of the system in 
the phase space are obtained. It is important to mention 
that if / is nondegenerate, the equation (|2.9[) has a very 
unique solution. The nondegeneracy of / means that the 
linear map t> : TM -> T*M defined by b{X) := b(X)f 
is an isomorphism. So, equation (|2.9[) can be solved 
uniquely for any Hamiltonian (X = \>~ l {dH)). On the 
other hand, the tensor has a zero-mode and a new con- 
straint arises, indicating that the iterative process goes 
on until the presymplectic matrix becomes nonsingular or 
singular. If this matrix is nonsingular, the Dirac brack- 
ets will be determined naturally. In 27(, the authors 
consider in detail the case when / is degenerated. 

The main idea of this embedding formalism is to intro- 
duce extra fields into the model in order to obstruct the 
solutions of the Hamiltonian equations of motion. We in- 
troduce two arbitrary functions that hinge on the original 
phase space and on the WZ variables, namely, ^(ai,Pi) 
and G(ai,Pi,rj), into the first-order one-form Lagrangian 
as follows 



V-(0) =y(°) + G( ai , Pi ,r,), 



(2.11) 



where the arbitrary function G(ai,pi, 77) is expressed as 
an expansion in terms of the WZ field, given by 

00 

G(a l ,p. l , V ) = J2G {n) (a l , Pt , V ),, (2.12) 

71=1 

where Q^ n \ca,Pi,Tj) ~ 77™, and satisfies the following 
boundary condition 



G{a l ,p i ,rj = 0) = 0. 



(2.13) 



The presymplectic variables were extended until they 
also encompass the WZ variable ^ e = (£ (0)e ,?7) (with 
6 = 1,2,..., 2N + 1) and the first-iterative presymplectic 
potential becomes 

00 

n=l 

In this context, the new canonical momenta are 



A (0) = j 4 0) > withg =1,2,..., 2N 
9 \ with 6= 2N + 1 

and the new presymplectic tensor is given by 



(2.15) 



m _ 



that is 



/AO) Aoy 

RO) _ [ Je/3 J6t] 



(2.16) 



(2.17) 



To sum up, basically, we have two steps: the first one 
is addressed to compute ^(ai,Pi) while the second one 
is dedicated to the calculation of G(ai,pi, 77). In order to 
begin with the first step, we impose that this new presym- 
plectic tensor (f^) has a zero-mode v, consequently, we 
obtain the following condition 



p(o)0/(°) = 



(2.18) 



At this point, / becomes degenerated and in consequence, 
we introduced an obstruction to solve the Hamiltonian 
equation of motion given by equation (|2.9[) . Assuming 
that the zero-mode v^ 6 is 



f,(°) = t.fi 



1) , 



(2.19) 



and using the relation given in (|2.18p combined with 
(|2.17[) . we have a system of equations, 



C^dt 



4 0) dC (0)e + ^dr ] -V ( ' 0) (0dt, (2-10) 



u e AO) , AO) _ 
M Je/3 + JriP — U ' 



(2.20) 
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where 



from 



AO) 



8Af 



(2.21) 



The matrix elements /i e are chosen in order to disclose 
the desired gauge symmetry. Note that in this formalism 
the zero-mode the gauge symmetry generator. It 

is worth to mention that this feature is important be- 
cause it opens up the possibility to disclose the desired 
hidden gauge symmetry from the noninvariant model. 
Prom relation (12.18[) some differential equations involv- 
ing ^{ai.pi) are obtained, i. e., the equation (|2.20[) . and 
after a straightforward computation, ^(ai,Pi) can be de- 
termined. 

In order to compute G(ai,pi,rf) following the second 
step of the method, it is mandatory that no constraints 
arise from the contraction of the zero-mode (pW> 8 ) with 
the gradient of the potential V^(ai,Pi,rj). This condi- 
tion generates a general system of differential equations, 
which is 



0. 



(2.22) 



dgW(ai,Pi,y) | 
drj 



f dgW(ai, P i,T)) 
d-q 



■ = 



(2.23) 



and that allows us to compute all the correction terms 
G^ n '{ai,Pi,Tj) as functions of 77. Notice that this poly- 
nomial expansion in terms of rj is equal to zero. Con- 
sequently, all the coefficients for each order in 77 must 
be identically zero. Given this, each correction term as 
function of 77 is determined. For a linear correction term, 
we have 



For a quadratic correction term, we have 







H(a i ,p i ,r))=vW(a i ,p i ) + G(a i ,p i ,r ] ), (2.27) 

and the zero-mode v^ e is identified as the generator of 
an infinitesimal gauge transformation, given by 



8~e = ev^ e , 



(2.28) 



where e is an infinitesimal parameter. 

In the following sections, we will apply the symplec- 
tic gauge-invariant formalism in some second-class con- 
strained Hamiltonian systems, Abelian and non-Abelian. 



III. THE ABELIAN PROCA MODEL 

Let us start with a simple Abelian case which is the 
Proca model whose dynamics is ruled by the Lagrangian 
density, 

C = - \f^F»» + l -m 2 A»A^ (3.1) 



where m is the mass, g^ v = diag(-\ ) and (from 

now on) F^y — d^Av — d^A^. As well known, the mass 
term breaks the gauge invariance of the usual Maxwell's 
theory. Hence, the Lagrangian density above represents 
a second-class system. 

To begin with the symplectic embedding procedure the 
Lagrangian density is reduced to its first-order form as 



£(°) = 7 r i i i -y(°), 



(3.2) 



where the presymplectic potential is 



- TO 2 Ai 

2 



1 



Ao(d i n i + ~m 2 Ao), 



W = + \f% 

(3.3) 

where 7Tj = Ai — diA$. From now on di — -^-r and the 
dot denote space and time derivatives, respectively. The 
presymplectic coordinates are ^i " 1 = (Ai, TTi, Ao) and the 
(2.24) corresponding one- form canonical momenta are given by 

0. The zeroth-iteration 



7T,- and a 



(0) 



(0) 



presymplectic matrix is 



) dg^(a i ,p i ,ri) , dgWfapuri) 



drj 



(2.25) 



From these both equations, a recursive equation for n > 2 
can be chosen as, 



1^ T" 



^(0)6 



0)] 



(2.26) 



which permits us to compute the remaining correction 
terms as functions of 77. This iterative process is suc- 
cessively repeated until (|2.22l) becomes identically zero. 
Consequently, the extra term G(ai,Pi,rj) is obtained ex- 
plicitly. Then, the gauge invariant Hamiltonian, iden- 
tified as being the presymplectic potential, is obtained 



/( o) = 



-Sij 



3' t 





I 6^(x-y), (3.4) 




which is a singular matrix. It has a zero-mode that gen- 
erates the constraint — diir l + m 2 A , identified as the 
Gauss law. We will include this constraint into the canon- 
ical part of the first-order Lagrangian in (|3.2|) intro- 
ducing a Lagrangian multiplier (/3). The first- iterated La- 
grangian can be written in terms of £^ = (Ai, m, Aq, /3) 
as = tt 2 Ai + Q/3 — with the following presym- 
plectic potential, 



V^ 1 



^ + Iff. + 1 m 2 ( Aq 2 + Ai 2) _ Aq ^ (3 5) 
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The first-iterated presymplectic matrix, computed as 



/ 





Sij 





\ 




Sji 








dV 













2 

m 


V 





-d x 


o 

— m z 


/ 



5^{x-y), (3.6) 



is a nonsingular matrix. Consequently, the Proca model 
is not a gauge invariant field theory. The Poisson brackets 
among the phase space fields can be obtained from the 
inverse of the presymplectic matrix. The Hamiltonian, 

n = vW\n=o = \nf - ^id%* j + \f% + \m 2 Al 

= \^My + lp* + (3.7) 

where the phase space metric is AO = g* — which 
completes the noninvariant analysis. 

The goal here is to disclose the gauge symmetry hid- 
den inside the model. Both arbitrary functions, and G 
rely on both the original phase space variables and the 
WZ variable (9). The former (^P) is introduced into the 
kinetic sector and the later (G), within the potential sec- 
tor of the first-order Lagrangian. The process starts with 
the computation of "J/ and it ends up the calculation of 
G. 

To clarify, the problem is that we have a second-class 
theory that is not (obviously) gauge-invariant. So, we 
have to make something if we want to transform this the- 
ory into a gauge invariant one. Thus, we have to intro- 
duce "something" new so that, interacting with the other 



fields and/or with itself, promotes a gauge-invariance 
"process" inside the theory. This "something" is the 
0-ficld and the "process" is the symplectic embedding 
fornalism. However, it is important to notice that this 
0-ficld, although now it is part of the symplectic coor- 
dinates of the extended phase space, it does not have a 
conjugated canonical momentum. Consequently, we be- 
lieve that, from the outset, this 9- field can not be inter- 
preted as the Stuckelberg field. Even if we consider that 
the canonical momentum conjugated to the Stuckelberg 
field is usually canceled by the Hamilton's equation of 
motion, the final result is different from the one obtained 
by the symplectic embedding. Note that the structures 
of the first-order Lagrangians for both methods are dif- 
ferent since in the Stuckelberg method [28| we have the 
pair lied and in the symplectic embedding we have ^/9, 
where iff is the arbitrary function introduced above. 

The first-order Lagrangian given in Eq. (13.21) . 

with the arbitrary terms, given by, 

£(°) =n i Ai+6(V + 7)-^, (3.8) 

where 

f(0) = h 2 + \ F * + \ m2 A < 2 + \ m2 Ao) 

+ G + ~fe 77 , (3.9) 

and * = ^(Ai,Tn,A ,9) and G = G{Ai,<Ki, A ,9) are 
the arbitrary functions to be determined as well as the 
constant k. Now, the presymplectic fields are — 
(Ai, iTi, A , 9, 7) while the presymplectic matrix is 



/ 




r(0) 



7TP 



V 



g 0% 8^{x-y) 





-<M (3) ( a 




dAf, 



y) 



OA- 

dAf 
9*3, 
dirf 



fS„9y 

6&(x-y) 







-S^(x-y) 




(3.10) 



with 



where # T 



je*e v g 



8%, 
89 v 



(3.11) 



0(x), 9 V = %), tt x = V(x) and % = *(y). 
Note that / is a 9 x 9 matrix with two space indexes in 
each entry. There is also an implicit time dependence, 
which comes from the coordinates and momenta. In the 
representation of /, described above, some zeros are ac- 
tually zero columns, zero lines or zero matrices. 

The corresponding zero-mode z^°)(x), the generator of 
the symmetry, satisfies the following relation, 



d 3 y^ } (y) f a p{x~y) = i) 



(3.12) 



The zero-mode does not generate a new constraint. How- 
ever, it determines the arbitrary function iff and conse- 
quently, it obtains the gauge invariant reformulation of 
the model. We will scrutinize the gauge symmetry re- 
lated to the following zero-mode, 



v= (1 lx3 lx3 b) 



(3.13) 



to be the zero-mode of j"' ', where b is a constant. 

The constraint generated by v is = —diir 1 — w? Aq. 
As we will see, v will produce a constraint which is equal 
to ft when 7 = 6> = 0. 

Following the procedure we know that v is a zero-mode 



Obtaining gauge invariant actions via symplectic embedding formalism 



7 



of /, one condition for \I/ is found, which is 
The constraint appears from the contraction, 



Now we add Xfl to C and consider A as a new indepen- 
dent field, that is, a Lagrange multiplier. Hence, 



(3.14) 



C {1) = it 1 At + (* + j)6 + XQ. — V. 



(3.16) 



d 6 y v a {x)^ 



SV(y) 



5£ a (x) 

-di-K -m A + / d y , A/ ^ 



(3.15) 



The presence of the constraint inside the kinetic part of 
the Lagrangian allows us to remove it from the potential 
part. Nevertheless, this common procedure will not help 
us here. Therefore no change happened in the potential. 

Setting ^ 1 > a = (A , A 1 , n l ,9, 7, A) as the new symplec- 
tic coordinates, where from now ona = 1,2, 10, and 
or, for short, fl = Q + Gq + bkj, where Gq is implicitly with the help of equation (|3.14|) . the symplectic matrix 



bkj, 



defined. 



is 



b5W 






9^ 

g£(g) 
8A3{y) 











' Sk-J (y) 



-MP) 

SA^x) 



Sn 1 (x) 



®xy 

5<® 











SA a (x) 



SGq(y) 



SGp(y) aVx(3) 



SG (y) 
88(x) 



bk5 (3 



SGo(x) 2t(3) SGo(x) c. x r(3) SG (x) SG (x) h i,:(3) 



(3.17) 



For the sake of clarity, it is convenient to use the notation 
5^ instead of 5(x — y). 

We can select two independent zero-modes to become 
the infinitesimal gauge generators, which are 

£(0) = (clq ad 1 cd l —kb l) , 

v M = (1 lx3 0i x3 b 0) = (D 0) . (3.18) 

The values of the constants ao, a and c can be freely cho- 
sen. We have to remember that different choices directly 
correspond to different gauge generators. As it will be 
shown, the value of b is also free. 

Naturally, other zero-mode structures are possible, 
some of which entail correspondence to both Wess- 
Zumino fields in each set. 

For i>( 7 ) just one condition is necessary to assure its 
zero-mode feature, namely, 



SG (y) 
5A (x) 



(m 2 -b 2 k)8(x-y). 



(3.19) 



This zero-mode is a generator of gauge symmetries. 
Therefore no new constraint may arise from its contrac- 
tion with the gradient of the potential. By equations 



(|3.16[) and (|3. 181) . we see that this condition is automat- 
ically fulfilled. 

There is a set of nontrivial equations that needs to be 
satisfied in order to V(q) be a zero-mode of Instead 
of evaluating them now, it seems to be easier to demand 
that may not give rise to a new constraint. Hence, 



0= / d 3 yi>? 6) (x) 



SV(y) 

= I d 3 y{a S{x~y){~d t 7r l -m 2 A Q ) 



(3.20) 



5A^(x) 



+ adl5{x-yWF t] -m 2 A t ) 
+ cd l J(x - y)(-n + d t A a ) + p 

+ Cd *8^{x) kb 86(x)J- 

The index x on d 1 means that the derivative must be 
evaluated with respect to x (i.e., d l x = d/dxi), and 

p^ = (a ,adi) (3.21) 



Equation (I3.20[) can be solved by considering G as a 
power series of 9 (and its spatial derivatives). Let Q n 
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be proportional to 9 n , so G — J2 n @n ■ The condition 
G(9 = 0) = leads to n > 1. Hence, 



Qi = — (-aofljTr* - m 2 p^A^ - cd% + cd^Ao). (3-22) 

The terms f?uf\ and ir^Sr do not contribute to the 
computation of Q\. But they do contribute to others 
GnS. After some straightforward calculations, one can 
find Q2 (without surface terms) as 



f 



2(kb) 2 



{c(2a + 0)8168*9 + m 2 /%e}. (3.23) 



The absence of and it 1 in Q2 implies that Q n = 
for all n > 3. Thus the function G is completely known 
and we can write down the expression for Go, which is, 

Go (2) = J d 3 y = ^(cd%6 - m 2 a 6). (3.24) 

Substituting this result in equation (|3.19j) . we have that 

2 

k = t^-. This fixes k as dependent of b. 

Our next and final step in order to prescribe the gauge 
embedded Lagrangian is to find ^. This can be done by 
demanding that urg-\ be a zero-mode of f^K Using (|3.24j) 
and (3.20) we have that, 



cd?5(x-y) + 



m 2 S^(x) 



b 5AJ(y) 



= 0, 



(3.25) 



ad*S(x-y) + 



m 2 Sty(x) 
b Stt^ (y) 



+ 8JS{x-y) = 0, (3.26) 



Also from the Euler-Lagrange equations, we have that 
'y = -^26. Thus, eliminating 7Tj and 7, the Lagrangian C 
can be expressed by 

~ I m 2 h ( 

C = - 7 F^F^ + —A»A^ + — \ -m 2 A 9 
4 2 m z I 

+ (1 - a)8 t 9{d % A -A l ) + aoSi&diAa - d,A l ) 

+ 6m 2 p^A t ,} (3.30) 



+ 



-{l-afdrOd 1 



8 t 98 1 9 



+ (1 - )(oo + c)a^a 4 + —p^op^t 



2m 2 



From the components of ve and i> 7 the infinitesimal 
gauge generators of the theory are obtained as 



SeA 

5 e A l 



sao — s, 
-ad l e, 



(3.31) 



where e is an infinitesimal time-dependent parameter. 

The symplectic formalism assures us that C is invariant 
under the above transformations for any constants b, ao 
and a (assuming they have proper dimensions, which are 
squared mass, mass and unit respectively). 

Usually, terms with more than two derivatives in the 
Lagrangian are not welcomed, these can be avoided by 
fixing o = l. 

If one wants an explicit Lorentz invariance, the con- 



stants need to be fixed as b = 



a = 1 and do = 



(alternatively, b could also be — m 2 ). With these values, 
the Lagrangian turn out to have a Stiickelberg aspect, 
that is 



SGpjx) 
59(y) 







(3.27) 



With equations p,14[) and (|3.25H3.27l) . up to an additive 
function just of 9 (action surface term), \1/ can be deter- 
mined as, 



r {m 2 A + cS i A i + (l-a)aV. i }. (3.28) 



Note that we can withdraw the term Af2 from £W 
without changing the dynamics. One can always redo 
the symplectic algorithm and find again the constraint Q. 
This will lead us back to C By varying C with respect 
to TTi and using Euler-Lagrange equations we find 

m = 8iA - At + -^{(1 - 0)8,9 + (oo + c)8 t 9}. (3.29) 

Note that the momenta are not the original ones (which 
are i^o), but when 9 is eliminated they are recovered. 



4 ^ 



-A^A^ - m 2 A^8^9 - 



-d»9d u 9. 
(3.32) 

The Lagrangian in Eq. (3.32) is not the most general 
one that can be achieved with the symplectic embedding 



method. Others structures of the zero-modes 



and 



&( 7 ) are also possible, and their components, together 
with the components of v 1 could also be field dependent. 

If we analyze this specific example (Abelian Proca 
model), we can see that the BFFT embedding in the 
Hamiltonian approach is the analogue of the Stiickelberg 
mechanism in the Lagrangian version. The extra field in 
the BFFT version is just the Stiickelberg field. This was 
discussed in (28|. This shows clearly what we said above 
about the fact that our 6>-field is not the Stiickelberg field. 

Following the symplectic embedding formalism, the 
zero-mode is the generator of the infinitesimal gauge 
transformation {SO = ev^), given by, 



SAfj, = -d^E, 



59 = -e . 



Indeed, for the above transformations the invariant 
Hamiltonian, identified as being the presymplectic po- 
tential V(°\ changes as STL = 0. It is a very easy task 
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to show the invariance of (3.35) under the above gauge 
transformations. Consequently, there is no need to carry 
out the Dirac analysis. 



IV. THE O(N) INVARIANT NONLINEAR 
SIGMA MODEL 

The O(N) nonlinear sigma model (NLSM) in two di- 
mensions is a free field theory for the multiplet o~ a = 
{a i, (72 1 • • • )ffn) satisfying a nonlinear constraint u\ = \. 
This model has its dynamics governed by the Lagrangian 
density 



1 



1 



£ = 2 d ^ ad ^a - 2 A (<7V„ - 1) 



(4.1) 



where // = 0, 1 and a is an index related to the O(N) 
symmetry group. 

The original second order Lagrangian in the velocity 
given in (|4.ip , is reduced into a first-order form, given by, 



£(0) 



(4.2) 



with = \ irl + \ A (o 2 a - 1) - \ a' a 2 . The presymplec- 
tic coordinates are t;^ = (o- a ,Tr a ,\). The presymplectic 



(4.3) 



tensor given by Eq. (|2.4p is 

/( o) = 



-Sab 

<5 ba | S(x-y). 




Since this matrix is singular, it has a zero-mode, = 
( 1 J . Contracting this zero-mode with the gradi- 
ent of the presymplectic potential , given above, the 
constraint obtained is, f2i = a\ — 1. The first-iteration 
Lagrangian is 



£« = n a a a + Qxp - V {1) |n 1= o, 



(4.4) 



with V 1 |n 1= o= 5 7Tq + |f'a- The presymplectic co- 
ordinates are ^i 1 ^ = (a a ,ir a ,p) with the following one- 
form canonical momenta, A^} — n a , A^} = and 



4 1 ' = K 



i 



The corresponding presymplectic tensor given by, 
= ( j 'M-r --//!. (4.5) 



-<5 a & CT a 

<5 ab | 

-<T 6 



is singular. Consequently, it has a zero-mode that gen- 
erates a new constraint, 0,2 = cr a ir a and the second- 
iteration Lagrangian is obtained as 



£(2) 



.p(al-l)+C(a a 7r a )-V^, (4.6) 



with F (2) = 

n 1= o- The enlarged presymplectic co- 

( 2^ 

ordinates are £q = (cr ,7r a ,p, C) and the new one-form 
canonical momenta are 

A%=n a , A%=0, AW=o*-l, A?=* a *». 



The corresponding matrix is 




<5(x - y), 



(4.7) 



which is a nonsingular matrix. This means that the 
NLSM is not a gauge invariant theory. Now, the original 
phase space will be extended with the introduction of a 
WZ field. Let us introduce ^(o- a , ir a , 9) and G(a a , Tr a ,0), 
into the first-order Lagrangian as follows, 



where the presymplectic potential is 



(4.8) 



1 



1 



+ £ A (o* - 1) + \ a' a 2 + G(a a ,ir a ,9), (4.9) 



with G((Ta, 7i"a, #) satisfying the relations given in Eqs. 
(12~T21 and (233) 1. 

The presymplectic coordinates are £a = (<7 a , 7r a , A, 0) 
with the following one-form canonical momenta, 



4 0) = 0. 
(4.10) 



The corresponding matrix /^°\ given by 



p(o) _ 



/ 





Sab 







£&a 



















V 


_a*x 








97r 6 





d*y \ 

da % 

Otv: 



fo 



6(x-y), (4.11) 



<} 



must be singular. This settles down the dependence re- 



lations between arbitrary function "J, namely, = 0, 
so, "J = ^(a a , Tr a) 6). This matrix has a zero-mode. 

Considering the symmetry generated by the following 
zero- mode, = (0a a 0ll Since this zero-mode 
and the presymplectic matrix (|4.11|) satisfy the relation 
(I3.12j) . ^ is determined as ^ — | a 2 a + c, where c is a 
constant parameter. 

We know that no more constraints are generated by 
the contraction of the zero-mode with the gradient of 
the potential. The correction terms as functions of 9 
can be explicitly computed. The first-order correction 
term in 9, Q^ x \ determined after an integration process, 
is Q^\a a , 7r Q , 9) = —a a ir a 6. Substituting this expression 
into Eq. (14.91) , the new Lagrangian is 



-1 



(4.12) 

However, the model is not yet gauge invariant because 
the contraction of the zero-mode with the gradient 
of the potential V produces a non zero value, indicating 
that it is necessary to compute the remaining correction 
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terms as functions of 9. It can be carried out just 
demanding that the zero-mode does not generate a new 
constraint. It allows us to determine the second order 
correction term given by = + ^a 2 9 2 . Substituting 
this result into the first-order Lagrangian (|4.12l) , we have 
that, 



£ (0) = vr„a fl 



2 a 



1 2 
- 7T„ 

2 



(4.13) 



Now the zero-mode does not produce new con- 
straints. Consequently, the model has a symmetry and 
all correction terms Q^ n ' with n > 3 are zero. 

We can recover the invariant second order Lagrangian 
from its first-order form. To this end, the canonical mo- 
menta must be eliminated from the Lagrangian (|4.13p . 
The canonical momenta are computed as 7r a = a a + o a 9. 
Inserting this result into the first-order gauge-invariant 
Lagrangian (|4.13[) . we have the second order Lagrangian 
as 

C =\d^a a - (J a a a )0 1)A, (4.14) 

with the following gauge invariant Hamiltonian, 

U=\rt + \a'l-{o a * a )6+\\(al-l) + \al9\ (4.15) 

From the Dirac point of view, VL\ arises as a secondary 
constraint from the temporal stability imposed on the 
primary constraints, ir\ and irg, and plays the role of the 
Gauss law, which generates the time independent gauge 
transformation. 

To proceed the quantization, we recognize the states 
of physical interest as those that are annihilated by fii. 

The infinitesimal gauge transformations 5£a = £v^\ 
are 



5o a = 0, 

5X = 0, 



5n a = so a , 

50 = E. 



Concerning these transformations the Hamiltonian 
changes as 5% = 0. 

Similar results were also obtained in the literature us- 
ing different methods based on Dirac's constraint frame- 
work dMi]- 

To disclose the hidden symmetry of the NLSM lying 
on the original phase space (<j a , 7r a ), we can use the Dirac 
method to obtain the set of constraints of the gauge 
invariant NLSM. It was described by the Lagrangian 
(|4.14p and Hamiltonian (|4.15[) . given by, 4>i = ir\, cf)2 = 
- 1), and <p\ = no and ip 2 = a a ir a -al0, where ir\ 
and ire are the canonical momenta conjugated to A and 
9, respectively. The corresponding Dirac matrix is sin- 
gular. However, there are nonvanishing Poisson brackets 
among some constraints, indicating that there are both 
second-class and first-class constraints. This problem is 
solved separating the second-class constraints from the 



first-class ones through constraint analysis. The set of 
first-class constraints is 



X2 



-(a 2 a -l) + n g , (4.16) 



while the set of second-class constraints is c^i = 7rg and 
?2 = CaTTa — <\fi ■ Since the second-class constraints are 
assumed to be equal to zero in a strong way [37j the Dirac 
brackets are constructed as 

{o-i(x),<7j(j/)}* = 0, 

{a i (x) ) n j (y)r = 5 ij 5(x-y) , (4.17) 
{n{x),irj(y)}* = 0. 

Hence, the gauge invariant Hamiltonian is rewritten as 



n = k 2 



1 (<Ja7T a ) 2 

2 a a a a 



2 a 



2 K 



(4.18) 



where the phase space metric My , given by My = <5y — 
^ 2 ^-, which is a singular matrix. The set of first-class 

constraints becomes 

1 



Xl = "A 



(4.19) 



Note that the constraint xi > originally a second-class con- 
straint, becomes the generator of gauge symmetries, sat- 
isfying the first-class property {\2,H} — 0. Due to this 
result, the infinitesimal gauge transformations are com- 
puted as 



5cr a 
5 ir a 
5X 



e{o- a ,X2} 
£{ir a ,X2} 
0. 



0, 



(4.20) 



where e is an infinitesimal parameter. It is easy to ver- 
ify that the Hamiltonian (|4.18p is invariant under these 
transformations because a a are eigenvectors of the phase 
space metric (My) with null eigenvalues. In this sec- 
tion we reproduced the results originally obtained in [38[ 
using an alternative point of view. 



V. THE GAUGE INVARIANT BOSONIZED 
CHIRAL SCHWINGER MODEL (CSM) 

It has been shown over the last decade that anoma- 
lous gauge theories in two dimensions can be consistently 
and unitarily quantized for both Abelian [j| |39|, H(| and 
non-Abelian [4l|, |42| cases. In this scenario, the two di- 
mensional model that has been extensively studied is the 
CSM. We start with the following Lagrangian density of 
the bosonized CSM with a > 1, 

C = ~ F^ + \ 8^d^ + q {sT ~ O d^A v 



(5.1) 
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,10 



Now, we have that g^ = diag(+l, —1) and e 
eio = 1. The Lagrangian in (|5.1I) . can be reduced into its 
first-order as follows, 

£<°> =7r^ + 7 r 1 i 1 -U<-°\ (5.2) 

where the zeroth-iterative presymplectic potential is 

+ q 2 A 1 + + qcj)'] 



Al [ - <77T0 



2 ? 2 (a+l)A x 



>'] . (5.3) 

The zeroth-iterative presymplectic variables are ^i ** = 
((f>, 7r^,, Aq, Ax, 7Tx) with the following one-form canonical 
momenta A a , 



A 



(0) 



4 



(0) 
Ai 



(5.4) 



The zeroth-iterative presymplectic tensor can be ob- 
tained as 



f {0) (x, y ) = 



/ -1 \ 
' 1 ' 

0-1 
\o 1 / 



S(x-y). (5.5) 



This matrix is obviously singular. Thus, it has a zero- 
mode that generates a constraint when contracted with 
the gradient of the potential U^°\ given by, 



fix = ^ 



dU (0) 



aei 0) 

= ir[ + q 2 (a - 1)A + q 2 A 1 + qn^ + qcj)', (5.6) 

that is identified as the Gauss law, which satisfies the 
Poisson algebra, {Vl-y(x), f2i(y)} = 0. The first-iterative 
Lagrangian £,W is 

=71-^ + 7^! + fi l9 ) -J7« , (5.7) 

with the first-order presymplectic potential given by 

- A [n[ + -q 2 (a - l)A Q + q 2 A t + qn ip + qcj)'] 



Ai[-qTt4- -q 2 (a + l)A 1 



(5.8) 



where EA 1 ) = U^°\ Therefore, the presymplectic vari- 
cj>, Aq,Ai,7Ti, r/) with the follow- 



(i) 



ables become 
ing one-form canonical momenta, 

A<P- 



n<t>, 



A 



4 l l - 4? = o, 



(5.9) 



AW = vri + q 2 (a - 1)A Q + q 2 A x + qic^ 






-1 











qdy 


\ 


1 














q 

















Q 
















-1 


q 2 













1 









qd x 


-q 


-Q 


-q 2 


-d x 





/ 



The corresponding matrix /W is then 
/ 

f (1) {x,y)= \[ \[ \[ n ^ S(x-y), 

(5.10) 

where Q = q"{a — 1) and this is a nonsingular matrix. 
This means that the model is not a gauge invariant the- 
ory. 

We have now two arbitrary functions 
\l>(0,7r0, A), Ai,7ri,0) and G(0,7r^,Ao,Ai,7ri,0), 
depending on both the original phase space variables 
and the WZ variable 9. The first-order Lagrangian can 
be rewritten as 



£ (0) = TcA + 7T\A\ + m - U (0) , 



(5.11) 



where 



1 



Ai K + 2 q2 ( a ~ + q2Al + ^ + q ^ 



Al [ - 
G(^,7T0, A , Ai,7Tl 



5? 2 (a + l)Ai 



(5.12) 



The enlarged presymplectic variables are now ^i '' = 
(</>, 7r^, Aq, Ai,tti, 9) with the following one-form canoni- 
cal momenta 



.4 



(0) 



and A 



.4 



(0) 



(0) 

Al 
* . 



7T1, 



.4 



(0) 



= = 



(5.13) 



The corresponding presymplectic matrix /(°) reads 



( 





-1 











a* H \ 




1 














di> y 



















a** 
a* B 
















-1 













1 













a* x 








\ 













S(x - y), 



where 



fo 



d9 x 

6{x), B y = 6{y), 



d9 v 



(5.14) 



(5.15) 



^(x) and * y = 

a** 



with 2 

Note that this matrix is singular since =0. So, 

* = *(0,7r^,Al,7Tl,0). 

We will now investigate the symmetry connected to the 
following zero-mode, P^ * 1 = (g — qd x 1 — q 2 —l), 
with bar representing a transpose matrix. 



Obtaining gauge invariant actions via symplectic embedding formalism 



12 



Now we have to multiply the zero-mode above by the 
presymplectic matrix (|5.14[) . Hence, some equations arise 
and after an integration "J is determined as \t = n[+q<j)'+ 
qn^ + q 2 A\. The presymplectic matrix (|5 . 14[) is rewritten 
as 



/CO) 



/ -1 qd v \ 

1 q 



-1 q 2 







ooi o a, 



6{x-y) (5.16) 



V -qd x -q -q 2 -d x / 



which is obviously singular and it has a zero- mode that, 
by construction, is given by the one above. 

The first-order correction term in0, £«,is determined 
by, 



6, 7T0, Al,7Tl, A 0: ( 



-^i6» + q 2 (a - 1)A[9 - q 2 6w 1 , 
(5.17) 



after an integration. Substituting this expression into the 
equation (|5.1ip , the new Lagrangian is 



with 



£ (0) = TT^ + TTlil + *<? - J7 (0) , 



Ao K + \q 2 {a ~ 1)^0 + q 2 M + qiTtp 
Axl-qnt- -g 2 (a + l)Ai-#'] 



(5.18) 



(5.19) 



The Lagrangian in f|5 . 18[) is not yet gauge invariant be- 
cause the zero-mode still generates new constraints, 
given by 



^^^(a-lf-^a-l)^^ 
The second order correction term Q^ 2 > is, 



(5.20) 



-q 2 (a-l)9 + q 2 (a-l)9" + g 4 



(5.21) 



50^ 
96* 



0. 



\ q 2 {a-l)9> 2 + \q±6 2 -\q 2 {a-l)e 2 



Hence, the first-order Lagrangian (|5 . 18[) becomes 

C {a) = tt^ + mAi + ^>6 - U m , (5.22) 



with the new presymplectic potential 

- A) [ttJ. + iq 2 (a - 1) A + q 2 Ai + 9 7r 



Ai [ - g7T0 - 
n x + g 2 (a 



- 1)6> A[ - q 2 9ir 1 



(5.23) 



1 



q 2 (a-l) 



1 



,4/)2 



-^ 2 (a-l) 



and the respective second-order Lagrangian is 



£ = 



1 



A' - A x 



0)0+ U + (Aq — Ax) 



Ai - 3Ai 



+ 9') 9 - 9' 9 - ~0' 2 
1 



(5.24) 



A [9-9')+aA 9+^(A 2 + A 2 ) 



+ LqZ(a-l)[Al-Al-6A' x ] . 

The contraction of the zero-mode with the new 
presymplectic potential above does not produce a new 
constraint. So, all correction terms with n > 3 are 
zero. The infinitesimal gauge transformations originated 
from the zero- mode (Sfy = ez/°)) are 



dip = qe, 
SA = e, 
Stti = —q e, 



Sn^ = qe', 
SAi = -e'. 
59 = -e. 



(5.25) 



It is easy to verify that the Hamiltonian, identified as 
being the new presymplectic potential U^ ', is invari- 
ant under these infinitesimal gauge transformation above, 
namely, SU^ = 0. 

We would like also to demonstrate that the anomaly 
was canceled. It will be carried out from Dirac's point 
of view. From the Lagrangian in Eq. (|5.22[) the sets 
of primary constraints are computed as, ifi — 7Tq, and 



-ne + ^ . The primary Hamiltonian is U 



(0) 



primary 



Xi 

U^ ' + \i<pi + JiXi- Since the constraint tpi has no time 
evolution, the following secondary constraint is calcu- 
lated as (fi2 = f^i — q 2 (o, — 1)9, and no more constraints 
arise from the temporal stability condition. In this way, 
the total Hamiltonian is 



U, 



(o) 

total 



U(°> + \i(pi + \ 2 (f2 + 7iXi- 



(5.26) 



The time stability condition for the constraint xi allows 
us to determine the Lagrange multiplier A3. In this way, 
the gauge invariant version of the model has three con- 
straints, namely, (pi, if 2 and xi- The corresponding Dirac 
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matrix, given by, 

/ -q 2 (a-l) 
C(x-y) = \ q 2 (a-l) q 2 (a~l) | S(x-y), 
\ -q 2 (a-l) 

(5.27) 

is singular. As the Dirac matrix is singular, the model 
has both first and second-class constraints. Through 
a constraint combination, we obtain a set of first-class 
constraints such as, xi = —ng + ty — iro, and a set of 
second-class constraints, given by (p\ — tpi and (f2 = 
ill — q 2 (a — 1)0- It is again easy to verify that \i 1S a 
first-class constraint, identified as the Gauss law, while 
the others are second-class constraints. Note that the 
anomaly was removed. Hence, the Gauss law is also rec- 
ognized as being the generator of the gauge transforma- 
tion given in Eq. (|5.25[) . 

The model has one first-class and two second-class 
constraints and the phase space dimensions result in 
eight dependent fields, i.e., ((f), 7r^,, Ax,ttx, Ao,Tro,9,ire)- 
The first-class constraint eliminates two fields, while the 
second-class constraints eliminate two fields. Hence, the 
model has four independent fields, i.e., there are two in- 
dependent degrees of freedom. 

In order to obtain the Dirac brackets, the set of second- 
class constraints, (px and </?2 , will be assumed equal to 
zero in a strong way. The Dirac brackets among the phase 
space fields are obtained as 



{<f>(x),w^(y)}* = 8(x - y), 

{<Kx),Mv)}*< ' 



-6(x - y), 



q(a-l) 

d y S(x-y), 



{ir^(x),A {y)}* 

{Ax(x),A ( y )y |, 
q 2 (a - 1) 

{Ax(x),itx(y)}* = S(x-y), 



q(a-l) 



{*i{x),Mv)V 



(o-l) 
{6(x),Tr e (y)}* = 5(x - y), 
{ir e (x), A (y)}* = 5(x - y), 



S(x - y), 



(5.28) 



the others are zero. Note that the Dirac brackets among 
the original phase space fields were obtained before. Af- 
ter this process, the model now have only one first- 
class constraint, x = Xi|yi=^ 2 =o = — + * identi- 
fied as the Gauss law, that satisfies the Poisson algebra, 
{x( x )i x(y)}* = 0- I n this way, the anomaly was elimi- 
nated, the symmetry is preserved, and the fundamental 
brackets among the original phase space fields were reob- 
taincd. Note that the Gauss law is the generator of the 
gauge symmetry given in (|5.25l) . 

Once more, the number of the independent degrees of 
freedom matches with the result obtained in the second- 
class case. The invariant model has a phase space 
((f), 7T0, Ax, 7Ti, 9, 7Tg), with six dependent fields, and has 



a first-class constraint which eliminates two fields. Con- 
sequently, the model has two independent degrees of free- 
dom. 

The remaining symmetry will be eliminated with the 
introduction of the unitary gauge-fixing term, given by 
6 = 0. Hence, both the noninvariant Hamiltonian and 
the corresponding Dirac brackets computed in the begin- 
ning of this section are reobtained, then recovering the 
anomaly. In this way, we conclude that the new symplec- 
tic gauge-invariant formalism does not change the physics 
contents inside the model. 



VI. THE NON-ABELIAN EXTENSION OF THE 
PROCA MODEL 

The non-Abelian extension of the Proca model has its 
dynamics governed by the following Lagrangian densitjf]], 



m 2 A a ^ a , 



(6.1) 



with F;„ = d„A a v - d v Al + gClA^Al. The antisym- 
metric tensor C£ c (C% c = ~C^ b ), are in fact a set of real 
constants, known as the structure constants of the gauge 
group, and satisfy the following property, 

The first-order form follows, 



£<°> = <i? - \(Kf + A a Q n a l -m 2 A^ a 



l m 2 A a A° a -^F^, 

-<6 



(6.3) 



where Q a = — gC b a nlA^ + m 2 A° a . The presymplectic 
variables are given by = (Af, irf, Ag) and the presym- 
plectic matrix is 

-5 ji 5 ba 0\ 

Q\8^\x-y). (6.4) 
n 0/ 



/ (0) = ( M ab 



Since this matrix is singular, it has a zero-mode that 
generates the constraint Q a - The first-order Lagrangian 
through a Lagrange multiplier, is, 



jrd) 



<i? + n a rf - i«) 2 + A a o n a - \m 2 A1A\ 



1 



^m 2 A a A° a 



(6.5) 



The new group of presymplectic variables is £° 
( Af , 7T? , Aq , rf ) , and the new presymplectic matrix is 

/ 



/ 



(i) 





















D ba (y) 











m 2 5 ab 


lf^(x) 


-D ba (x) 


-m 2 5 ba 






S^{x-y) 



(6.6) 



1 The BFFT embedding for the non-Abelian Proca theory was 
discussed in [36(| 
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where D ab {w) = 8 ab df - gCfAf(w). This matrix is 
nonsingular. The model now can be reformulated as a 
gauge invariant field theory. 

The first-order Lagrangian (|6.3|) can be rewritten as 



£(0) 



■k 1 A a 



(6.7) 



G, 



where the arbitrary functions are 
* a = * a (A?,Tr?,A%,9 a ), 



(6.8) 



G ee G(Af, *?, Ag, 9 a ) = J2 G n (Af, n?,A%, 9 a ) . 

n=0 

The G function obeys a boundary condition given by, 

G ee (A?,ir?,A$,9 a = 0) = g°(At,^,A a ,e a - 0) = 0. 

(6.9) 

In this context, the corresponding presymplectic matrix 



/( o ) = 



-5 H 5 ba 



S l3 5 ab 




— m \ 



J 



dA°;(x) 
8A%(x) 

d<H a (x) d^g(x) d^q(x) n 
\ dA^y) drfiy) dAl(y) 

(6.10) 

In order to determine the functions ~^! a , we ana- 
lyze the symmetry related to the following zero-mode, 
= (df l) which produces a set of differential 
equations which allows us to compute the ^ a function 
as ^ a = —diTT l a {x). Consequently, the first-order La- 
grangian can be rewritten as 



£(°)=<i»-(a I <)0»^(°) 

where the presymplectic potential is 

^ (0) = \{<f - A a n a + \ m 2 A1A\ + \ 
1 



(6.11) 



™ 2 A«Al 



G. (6.12) 
For the hidden symmetry inside the model, we have, 



The linear correction term in 9 is given by, 



(6.13) 



df 



= 0. 



dA{(w) d0f(w) ^ 
After an integration we have that 

g (1 \x) 

= - 9 C h fa dUA a (x)4(x))d f (x) - m\d*A))9t\x) 



(6.14) 



dF^jx) 
dA{(y) _ 



9f(y). 



3.15) 



Now, we will compute the quadratic term, namely, 

,dg^(x) gW(x)~ 



8A{(w) d0f{w) 



0. 



(6.16) 



Integrating this relation in 9f(w), the quadratic correc- 
tion term is obtained as 



gv\x) = l -m 2 {dief{x)Y 



0f(x) Jw 



[(df^} b )e"(y)] , 



(6.17) 



where 



« _ dF Sj {x)dF«(x) p . ■ „ 

}b dA{(w) dA\{y) k > y dA{(w)dA b (y) 

(6.18) 

In this way, two correction terms as functions of 9 a 
(g^ (x) and g^ (x)) were computed yet. Let us compute 
the first one. It can be carried out from the following 
relation, 



g^(x) 



d9 g {z) 
and we can write that, 

g (3) (x) 

' ' ' di 



0f(x) Jw 





d? / d, 



dA k f 



kl 

fb a b 



dA 9 n (z) 



(.'/) 



(6.19) 



0a(z) J z 



s (x) Jw 



Qjkl 

fl? / df^^9 b (y). 



dA?(z) 



Finally, the last correction term is, 



f f dl 


I I 9; -J 


I0 h (v) Jv 


j0g(z)Jz J0 



B w / gy °_ A fb 

k J y 1 dAi(v)dA-(z) 



\y)- 



(6.20) 



Therefore, the gauge invariant first-order Lagrangian is 

£(°)=<i?-(5 i <)^-F(°), (6.21) 

where the presymplectic potential, identified as being the 
gauge invariant Hamiltonian, is given by 

VM = ~(ni) 2 -A a n a + ~m 2 (A?) 2 



+ \m 2 {A a Q ) 2 + \{F% j ) 2 

- gC» fa df (A a (x)4(x))9f(x) 

- m 2 {&?A\)9 f {x) + -m 2 {d l 9 f {x)) 

J 2 

+ h+h+h + h 



3.22) 
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where 



dF*1{ 3 



e f (y) 



kj ^ J dA{(y) j 
h = +\e f {w) I df j [(d?A%)9 b (y)} 



h = - 



h 



9s(z) 

1 

2 



0,. 



f(x) Jw 



1 S d * 


I J 9 *J 


>e h {v) Jv 


J9 g {z)Jz Je 



f)A kl 
v 1 dA-(z) 



\y) 



df- 



d 2 Af b 



,„ dAKv)dA^(z) 
The second-order Lagrangian is 

2 



m = \ (At 



2A a l d l 9 a + Afd l A 0a + d t 9 a 



- \ (d t A a f - AfdX - 8^8% 

+ gCt c d l A 0a A b Ac - ZgC^A^Al - gC a hc A l a A b Al 

- gC% c d%A\Al + \g 2 {C£ c A b A c ) 2 - ^g 2 (C&flj^) 5 
-gCZ c A c d l 9 b (A* 
+ 1 -m\Alf -\m\Atf 



+ diOa + g*CS c C amn a i ff'AZA$A m 



I ( F a ) 2 
4 v v> 



(6.23) 



'i 2 (diAl) 



x - m 2 (d l e a ) 2 



h 



Let us make an analysis using the Dirac point of 
view. We start with the set of primary constraints, 
f2° = d l irf + TTg and Xi — ^o- For ^ ne nrs t se * °f con ~ 
straints, the time stability condition is satisfied (fif = 0). 
For the second one, the following secondary constraints 
are required, \2 — ^ a ~ 9C ba "^l^i^ ■ Hence, the total 
Hamiltonian is U = U + A*fiJ + C^Xi + Cxh wh ere A*, 
and Ca are Lagrange multipliers. Since the Poisson 
brackets among those constraints are 

{£!?(»), nJCtf)} = {n$(x), x b 1 (y)} = o, 

{nUx), X b 2 (y)} = 0, {rf(aO,X2(iO} = -mW 3 >(a?-0), 

{xSO^xSfo)} = 2 5 C7f X ^(x)<5( 3 )(f-y) 

-2p 2 CM(^ (3) (^9), 

no more constraints arise. Notice that some brack- 
ets above are zero, indicating that there are both first 
and second-class constraints. Indeed, the first-class con- 
straint is £2" and the second-class are Xi and xt- ^ n 
agreement with Dirac's procedure, the second-class con- 
straints can be taken equal to zero in a strong way. This 



allows us to compute the primary Dirac brackets. Due to 
the Maskawa-Nakajima theorem [37], the primary Dirac 
brackets among the phase space fields are canonical. To 
demonstrate this, the brackets are computed explicitly. 
The Dirac matrix is 



C = 





m 2 S dc 



-m 2 6 cd 

B cd 



(6.24) 



with B cd = 2gGl d x\{x) - 2gm 2 C£ d A b (x). The inverse of 
the Dirac matrix is 



C (-D 



1 



B'' 



6 cd 




6&(x-y). 



(6.25) 



In accordance with Dirac's process, the Dirac brackets 
among the phase space fields are obtained as 



{At(x),n b (y)Y 
{A1{x),A b {y)Y 



S ab S^\x-y), 
-l-dfS^(x-y) 



l^CfA{(x)5^(x-y), 



{n?(x),A b (y)}* = -±^gC?nf6&(x-y), 
{A a (x),A b (y)r = -^CfA e (x)SW(x-ff), 



{A a {x)y e {y)Y 



,Cfdfn?x)5^(x-y),(^) 



{9 a (x),4(y)r =S ab S^(x-y), 

and the others are zero. Finally, the infinitesimal gauge 
transformations are obtained, namely, 



5 At = ~d*e a , 
Snf = 0, 

6A% = 0, (6.27) 

50 a = e a , 

K = °. 

which lead us to the invariant Hamiltonian. 

To demonstrate that the gauge invariant formulation 
of the non-Abelian Proca model is dynamically equiv- 
alent to the original noninvariant model, the symme- 
try is fixed by using the unitary gauge fixing proce- 
dure, ip a — 9 a w 0, which leads to the bracket below, 
{toi(x),<p b (y)} = ~5 ab S^(x-y). So, a new Dirac brack- 
ets must be computed. The corresponding Dirac matrix 
for this set of constraints is 



C = 



-1 

1 



5 (3) (x~y). 



(6.28) 



Using the inverse of this matrix, the Dirac brackets 
among the physical phase space fields can be computed, 
which is equal to the one calculated from the original 
description. This result demonstrate again that the sym- 
plectic formalism can be seen in fact as a mapping be- 
tween the original theory and the final one. 
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VII. HIDDEN SYMMETRIES OF A FLUID 
DYNAMICAL MODEL 

It was demonstrated in [43( that the relativistic the- 
ories of membranes are integrable systems through the 
transformation of the problem into a two-dimensional 
fluid dynamics one. In this case, the potential term is 
proportional to the inverse of mass density, i.e., V oc p^ 1 . 
This subject is connected in some way to other fields 
like the parton model 14411 , hydrodynamical description of 
quantum mechanics [451 ] . black hole cosmology [46j and 
hydrodynamics of superfluid systems (47| . 

The majority of these cases aimed to find the solutions 
of the Galileo invariant system in d-dimcnsions in connec- 
tion with the solutions of the relativistic d-brane system 
in (d+ l)-dimension (48j . 

Let us start with the linear (or nonlinear) , Schrodinger 
Lagrangian theory defined in a d-dimcnsional (r) space 
evolving in time (t), 

L s = [d d r[W*V- i (W)-(W) - V(***)] (7.1) 



where V represents any nonlinear interaction and 

* = p 1/2 e lB , (7.2) 



where p = p(t, r) and 9 = 8(t, r) is the velocity potential 
[47] | into the Schrodinger Lagrangian. Now we can write 
the fluid dynamical model 



L = J d d r[6 p - ipV0 -V0 - V{p)} 



(7.3) 



where V{p) — V(p) + 



1 (vpr 

p 



which is the hydrody- 



namical form of the Schrodinger theory (details can be 
found in |50|). We can say here that the connection be- 
tween the fluid model, the membrane and its general- 
ization to d-brane systems only appears under the very 
specific density-dependent interaction potential V = g/p 
[481 ] . The symmetries and the corresponding generators 
for the fluid dynamics model can be found in [50|. 

The Lagrangian in (I7.3[) is already written in first- 
order form, so we can write C(°> — Op— where 
= \pdi6d l 9 + V(p). The symplectic coordinates 



(0) 



9 and A 



(0) 



0. 



The zeroth-iterative symplectic matrix 



f 



(0) 



5(d) ( r 







- < y(«0(r-r / ) 




(7.4) 



is non-singular and the model is not gauge invariant. 
We have that, 



£ (0) = 0p + - V {a) 



(7.5) 



where 



\pd l 9d l 9 + V(p), and # = V(p,6) and 
G = G{p, 9, rj) and the symplectic matrix is 



f(0) 



/ 

5* r 



SpJF) 





'Se(r') 



8p(r) 

8 <b T , 
SWT) 




(7.6) 



where * r = *(p(r)), 9(r)), * r , = V(p(r'),0(r')) and the 
zero-mode (r) satisfies the relation 



d d vv^{v)f~ e - {v,v') 



(7.7) 



which provides us with the set of equations that allows 
the determination of \& explicitly. 

The WZ gauge symmetry is related to the following 
zero-mode = (11 — 1) since this zero-mode and the 
symplectic matrix (|7.6[) must satisfy the gauge symmetry 
condition in (|7.7[) . We can obtain the set of equations 



fly-Uid) ( r -r') 



5^ r 







^ r (_ jM(r _ r0 + ^_) = „ (7 . 



.6p{r') + 6p(r') 



= 



After an integration, \1/ can be written as ^(r) = 9(r) 
p(r) and the symplectic matrix, 

-1 -1 

/W =(l0 1 ) <5 (d) (r - r') (7.9) 

1 -1 

which is singular. Hence, the first-order Lagrangian is 
= 9p + (9 - p)f] - V<® , (7.10) 

where V m = \ (p + rf) {d r 9f + V(p - 77). The zero- 
mode v 1 ^ is the generator of infinitesimal gauge trans- 
formation, then 

5p(r,t) =e(r',t)5 ( - d) (r-r') 

59(r,t) =e(r',t)S {d) (r-r') (7.11) 

Sri(r,t) = -e(r',t)S (d) (r-r') 

where s(r,t) is an infinitesimal time-dependent parame- 
ter. It was shown in [50( that the Lagrangian density in 
([7.10P becomes, 

C = -(p-v)0-l(p-v)(dief +V(p-ri) , (7.12) 

which is the same result obtained in [49| . This La- 
grangian can also be written as 

£=(P-V) dViP g ~ V) +V(p- V ) or 

C = -(p-r,) dViP d ~ r,) +V(p-r,) , (7.13) 

which depends on both equations of motion for rj and for 
p respectively. 

As a second example [5(| of zero-mode for this theory 
we have that the zero-mode = (1 — 1) will generate 
the gauge-invariant first-order Lagrangian [50l |. 

I = -(p + n) 9 - \ ( P + v ) (di 9f + v( P + n) , (7.14) 
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and the infinitesimal gauge transformations will be ob- 
tained as 

Mr.*) =e(r',i)<S (d) (r-r') , 56(r,t) =0 

6 rj(r, t) = -e(r', t) 5 {d) (r - r') (7.15) 

In |50( | , the interested reader can obtain other six sym- 
metries within the fluid model described in (|7.ip . 



VIII. HIDDEN SYMMETRY IN THE 
ROTATIONAL FLUID MODEL 

In some of us have demonstrated that the irrota- 
tional fluid model has a set of dynamically equivalent WZ 
gauge invariant versions. Also, the extra global symme- 
tries, namely, Galileo antiboost and time rescaling, first 
obtained in [5l|, were promoted to local symmetries. 

In this section, we will investigate the symmetries of 
the rotational fluid model^ but now they have an extra 
term, like kp(di9 + adi/3) , where k is a constant, p = 
p(t,f) is the mass density and 9 = 9(t,f) is the velocity 
potential. The terms introduce a dissipative force into 
the model. 

To establish our ideas, we will investigate how the in- 
clusion of a dissipative term affects the dynamics of the 
fluid model considering only one of these zero modes in- 
troduced in |50| |. 

It is well known that systems that have vorticity 
and/or viscosity have Casimir invariants which obstruct 
the construction of a canonical formalism for a fluid, as 
demonstrated in 



52|. 



However, this obstruction can be 
eliminated using the Clebsch parameters, as it was shown 
by Lin [53[ and by two of us in 54 1. In fact, with the 
introduction of Clebsch parameters, it is possible to ob- 
tain a Lagrangian density for the rotational fluid with 
dissipation, 3-dimensional, as being 

C = —p{9 + a/3) - V, (8.1) 

where the presymplectic potential is 

V = -(1 - k)p{d t 9 + ad l /3)(d l 9 + a&p) + V{p). (8.2) 

The presymplectic coordinates are = (p,6,a,j3) 

and the corresponding zeroth-iterative one-form canoni- 
cal momenta is given by A p ^ = 0, — —p, Aa^ = 
0, A^ = —ap. The zeroth-iteration presymplectic ma- 




-S(f-r*) -a 

5(r-r*) 

-p 

a p 



(8.3) 



is a nonsingular matrix. The model is not gauge- 
invariant. 

Although the symplectic formalism does not restrain 
the dimension of the model, we choose a 3-dimensional 



description for the rotational fluid to place our work in a 
correct perspective in comparison with others. 

The first-order Lagrangian £(°>, equation (|8.1|) . with 
additional arbitrary terms (vP, G) is given by 



£ (0) = -p(9 + ap) + 4>r/ - V (0) , 



(8.4) 



where 

V (0) = ~(1 - k)p{d,9 + ad^)(d l 9 + ad l (3) + V{p) + G, 

(8.5) 

and ^ = ^(p, 9) and G = G(p,9,r]) are arbitrary func- 
tions that must be calculated. The presymplectic coor- 
dinates are |w = (p,6,a, /3,r]) while the presymplectic 
matrix is 



?(0) 






-Sir-r 1 ) 





—a 


S(r-7) 




















-P 








P 













Sp(r') 


S9(r>) 


Sa(f') 





5p(r, 
S0(r) 
6a(r) 





where = ^(p(r),9(r)) and $p = *(p(r*), 9(7*)). 

With a general zero-mode, = ( a b c d — 1 ), we 
have the following set of differential equations 



dr 
dr 
dr 
dr 



bS(f- f) + ^+ daS(f- r ) ) = 0, 
op 1 



—a5(f— r ) + 



0. 



(8.6) 



2t+dp5(r-7?) =0, 

^ - aa - cpS(f- fu = 0. 



Hence we obtain that 

\&(r) = — bp — adp + a9 + aa/3 + cpfi. 



(8.7) 



In order to have a final solution, we consider that a = c = 
0, then i/W = ( 1 1 -1 ), and *(r) = -(1 + a)p. 
Then, 



E 



dr v K 

dr' {(1 + a)p{d[9 + ad'^d'^f - f) 

( dQ {n) dg^ dg^ n) ' 



V 89 



<9/3 



drj 



")} 







(8.8) 



To compute the first correction term as function of 77, 
0^, we pick up the terms in equation (|8.8[) with zeroth- 
order in rj, thus 



df 1 



{(1 + a)p(die + adim^ - f) - 



drj 



= 



(8.9) 
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where d[ = ^ • The linear correction term as function of 
X] is = (1 + a)p (di9 + adi/3) dif]. For the quadratic 
correction terms in equation (|8.8p . we have that 



dgw dgw egw 



89 



d/3 



0. 



(8.10) 



And the second-order correction term is — 
— — —pdi'qd 1 ?]. For the cubic correction terms in equa- 
we can write 



2 

tion 



dgW ogm d g(3) 



m 



9/3 



dr\ 

dg^ 

drj 



= 0, 
= 0, 



3.11) 



which allows us to conclude that g^ = for n > 3. 
Hence, the gauge-invariant first-order Lagrangian is writ- 
ten as 

£W = -p{9 + a$) - (1 + a)pf] - V {0 \ 

= -p(6 + r))-ap(l3 + r))-VW (8.12) 

where the presymplectic potential is 

V {0) = -(1 - k)p{di9 + ad l /3)(d i 9 + ad l /3) 
+ (1 +a)p(di9 + adi/3) diV 

+ ^¥-pd i r,d i r, + V(p) ) 

= ±(l-k)p[d i (e + r,)+ad i (f3 + r l )] 

x [d l {9 + tj) + ad l {p + r})] + V(p) . (8.13) 

The Lagrangian (|8.12[) can also be written as C^> = 
p after solving the equation of motion for p. 

The zero-mode v^ 1 is the generator of infinitesimal 
gauge transformations (SO = ei>^). Then, 

5p{r.t) = 0, 

59(r,t) =s(r,t) £(f - f*), 
8a(r,t)=0, (8.14) 
5/3{f,t) =e(r,t) S^-r 1 ), 
5i](r, t) = — e(r, t) S(r — f 1 ), 

where e(r, t) is an infinitesimal time-dependent param- 
eter. In fact, under the infinitesimal transformations 
above, the invariant Hamiltonian (V^) changes as 
SV^ = 0. Considering the following transformations, 
9 = 9 + r) and (3 — f3 + 77, then the Lagrangian density, 
equation (|8.12[) . and the Hamiltonian, equation (|8.13p . 
become 



£(°) = -p(9 + ap)-V^\ 

V {0) = \ P {di9 + adiPWO + ad l ~P) + V(p) 



.15) 



These expressions are identical to the original expressions 
for the Lagrangian in equation (|8.1|) and the Hamiltonian 
in equation (|8.2j) . respectively. 

Thus, at this point, it is important to point out that 
exist a hidden symmetry into the rotational fluid model. 
For the set of differential equations obtained we have no 
other solutions, i.e., the model has, in fact, only one hid- 
den symmetry. Based on the investigation done by some 
of us in Ref . |50l | . where the extra global symmetries 
proposed in [5l| are promoted to local, we can conclude 
that these extra global symmetries do not exist in the 
rotational fluid model. 



IX. FINAL DISCUSSIONS 

From the Dirac point of view, a system classified as 
a gauge invariant theory is one that has first-class con- 
straints. When a theory has second-class constraints, 
the gauge invariance can be recovered by converting the 
second-class constraints into first-class constraints. In 
the literature there is a great variety of techniques, with 
pros and cons, to promote this kind of conversion. 

In this work we are concerned not only with the gauge 
invariance of second-class systems but also with the ob- 
tainment of a theory physically equivalent to the original 
one. We will see why we believe that symplectic em- 
bedding formalism is the most adequate technique [T2I ]. 
Besides, as demonstrated by some of us, this method has 
the advantage that a convenient choice of a convenient 
zero-mode can lead to a theory physically equivalent to 
the original one through the elimination of the WZ terms 
0, H[. This means a new interpretation of the method. 
However, as recently demonstrated by some of us [56| the 
choice of the zero-mode must obey some "boundary con- 
ditions" . In other words, we can say that the physical 
coherence must guide us to choose the correct (or con- 
venient) zero-mode. But at the same time it is possible 
to obtain, as mentioned just above, a whole family of 
physically equivalent actions. 

Speaking in another way, we can say that this par- 
ticular mapping between the original action and the re- 
spective final gauge invariant theory can be interpreted 
as a kind of physical equivalence. Whenever necessary 
or convenient, the gauge invariance of the final actions 
obtained here was demonstrated via the Dirac analysis. 
The characterization of the final action as a first-class 
system corroborates the success of this process of sym- 
plectic embedding. 

Firstly in this paper we used a kind of toy model, the 
Proca model, to illustrate the procedure. After that, we 
apply the formalism to the non-linear sigma model and 
to the chiral Schwinger model. In the NLSM, a hidden 
symmetry lying on the original phase s pac e was disclosed, 
differently from other approaches (l2l. |30L l3~i| , where the 
symmetry resides on the extended WZ phase space. 

In the CSM, the chiral anomaly was eliminated and 
the gauge symmetry was recovered. It is important to 
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notice that this result was achieved introducing one WZ 
field while other schemes in the literature thrive with the 
introduction of two or more WZ fields, which is the origin 
of the ambiguity problem. 

Besides, we showed in the context of a non-Abelian 
model (the non-Abelian Proca model) that the symplec- 
tic embedding formalism can be used without any restric- 
tions concerning the noninvariant model algebra. Other 
constraint conversion techniques work since the algebra 
was previously and necessarily taken into account. 

We also have brought a gauge-invariant version for 
the rotational fluid model. As a consequence, we have 
demonstrated that the hidden symmetry found is unique. 
Although we have studied the rotational fluid model 
with an extra term, which introduces dissipation into 
the model, the results are also valid without viscosity 



(k = 0). We have noted that, although we have dissipa- 
tion, this fluid model has a hidden symmetry, which does 
not belong to the other group of symmetries obtained for 
the irrotational fluid model [50fl. So, the local version 
of the extra global symmetries [51j does not exist in the 
rotational fluid model, with dissipation or not. Further- 
more, the physical meaning of the hidden symmetry can 
be interpreted. 
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